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Abstract 

The quantum interference effect on the quasiparticle density of 
states (DOS) is studied with the diagrammatic technique in two-dimensional 
(i-wave superconductors with dilute nonmagnetic impurities both near 
the Born and near the unitary limits. We derive in details the expres- 
sions of the Goldstone modes (cooperon and diffuson) for quasiparticle 
diffusion. The DOS for generic Fermi surfaces is shown to be subject 
to a quantum interference correction of logarithmic suppression, but 
with various renormalization factors for the Born and unitary limits. 
Upon approaching the combined limit of unitarity and nested Fermi 
surface, the DOS correction is found to become a ^-function of the 
energy, which can be used to account for the resonant peak found by 
the numerical studies. 

PACS numbers: 74.25.Jb, 71.20.-b, 73.20.Fz 

I. Introduction 

Since the experiments revealed the d-wave symmetry of the order parameter in 
cuprate superconductors the physics of low-energy quasiparticle excitations in dis- 
ordered two-dimensional (2D) d-wave superconductors has been a subject of ongoing 
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intensive research |2). The characteristic feature of the d x 2_ y 2-wave pairing state is the 
existence of four nodal points where the order parameter vanishes. In the vicinity of the 
gap nodes there exist low-lying Dirac-type quasiparticle excitations. An understanding 
of the disorder effect induced by randomly-distributed impurities on these low-energy 
states is essential for the elucidation of the thermodynamic and transport properties 
of disordered d-wave superconductors. During the years a number of theoretical ap- 
proaches, such as the self-consistent approximation schemes [3-9], non-perturbative 
methods [10-17], and numerical studies [18-22], have been developed to calculate the 
quasiparticle density of states (DOS) in the presence of disorder. Unfortunately, many 
of these theories contradict each other. The DOS at zero energy was variously pre- 
dicted to be finite [3-10], vanishing [11-14,18,21], and divergent [15,20]. Recently, it 
has been made clear that [16,17,19,22] much of this controversy roots in the fact that 
the d-wave superconductor is fundamentally sensitive to the details of disorder, as well 
as to certain symmetries of the normal state Hamiltonian. 

In order to clarify the physics of the various asymptotic results for the DOS, 



Yashenkin et al. [23] analyzed the diffusive modes in disordered d-wave superconduc- 
tors and calculated the weak-localization correction to the DOS with the diagrammatic 
technique. It is well known that the weak localization in electron systems is a man- 
ifestation of the quantum interference (QI) effect, which stems from the existence of 
the Goldstone modes (cooperon and diffuson) p3|. As pointed out by Altland and 



Zirnbauer [25|, the Goldstone modes in a superconductor have different features from 
those in a normal metal, for the local particle-hole symmetry of the superconducting 
state gives rise to a combination of impurity- and Andreev-scattering processes. As 
a result, every cooperon or diffuson mode in the retarded-advanced (RA) channel en- 
tails a corresponding mode in the retarded-retarded (RR) or advanced- advanced (AA) 
channel. In the combined limit of the unitarity and nested Fermi surface (the UN 
limit), each of these 0-mode cooperon and diffuson has a 7r-mode counterpart due to 
the global particle-hole symmetry p3|, p6|]. For disordered d-wave superconductors, it 
was found p3] that in generic situations the quasiparticle DOS is subject to a loga- 
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rithmic weak-localization correction due to the 0-mode cooperon, and the existence of 
diffusive 7r-modes can give rise to a finite, or even divergent zero-energy DOS upon 
approaching the UN limit. 

The weak-localization calculations carried out by Yashenkin et al. appear to cap- 
ture the physical origin of the discrepancies between predicted low-energy quasiparticle 
DOS. Furthermore, the quasiparticle weak-localization effect was also suggested to have 
important influences on transport properties such as the electrical Q, spin ||13|| , ther- 
mal and optical 2£] conductivities of d- wave superconductors. While the QI effects 



have been widely investigated for disordered normal metals [24|, a corresponding the- 
ory for random Dirac fermions in superconducting cuprates is far from well established, 
and thus highly deserves further development. 

This paper presents an intensive study of the QI effect on the quasiparticle DOS 
in 2D (i-wave superconductors with dilute nonmagnetic impurities both near the Born 
and near the unitary limits. These two limiting cases are considered to be closely 
related with the disorder effects in cuprate superconductors. It is reasonable that 
the disorder due to defects off the copper-oxygen plane may be treated in the Born 
approximation and that defects in the plane may be in the unitary limit ||. Albeit 
sharing certain aspects with Ref. p3[ , we further develop the weak-localization theory 
in d-wave superconductors, and obtain some new results in this paper. First, the 
expressions of the Goldstone modes for quasiparticle diffusion are derived in details 
both near the Born and near the unitary limits. Second, we calculate the additional 
contributions to the DOS from those lowest-order self-energy diagrams containing non- 
singular ladders, which were not taking into account previously. These diagrams are 
shown to give rise to various renormalization factors for the DOS correction in the 
Born and unitary limits. Third, by taking into account a new nontrivial self-energy 
diagram with the 7r-mode diffuson, we show that the QI correction to the DOS becomes 
a (^-function of the energy upon approaching the UN limit. This result can be used to 
account for the resonant peak found by the previous numerical studies Jl9| , p0fl . 

The structure of this paper is as follows. In Sec. II, the commonly used self- 
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consistent T-matrix approximation (SCTMA) is described for a weakly-disordered 
d x 2_ y 2-wave superconductor. Using the SCTMA, we derive in Sec. Ill the expressions 
of the 0-mode and 7r-mode cooperons and diffusons. The QI correction to the quasipar- 
ticle DOS is calculated in Sec. IV, and the conclusions are summarized in Sec. V. The 
Appendix provides some mathematical formulas, which are useful for our derivations. 



II. SCTMA for d-wave superconductors 

Let us consider a most extensively studied model for a 2D d x 2_ y 2-wave supercon- 
ductor. In the tight-binding approximation, the normal-state dispersion of a square 
lattice is given by = —t(cosk x a + cos k y a) — [i where o is the lattice constant, t 
is the nearest-neighbor hopping integral, and /j, is the chemical potential. The nested 
Fermi surface corresponds to the half- filling case (// = 0). The order parameter of 
the d x 2_ y 2-w&ve pairing state can be expressed by A& = Ao(cos k x a — cos k y a). The 
four gap nodes are given by k n = ±(ko,±ko) with ko = iX/a) arccos (— fJ>/2t), which 
satisfy £fc n = Afc n = 0. In the vicinity of these nodes, the quasiparticle spectrum 



can be linearized as ej. = J + A| ss y (vfk) 2 + (v g -k) 2 , where Vf = (d£ k /dk) kn , 
v g = (dA k /dk) kn , and k is the momentum measured from the node k n . A direct 
calculation yields vj = v g 

In the presence of randomly-distributed nonmagnetic impurities, the time-reversal 
and spin-rotational symmetries remain preserved. Thus the system belongs to symme- 
try class CI in the classification of Ref. [25|. The impurity potential is assumed to be 



point-like, meaning that the intra- and inter-node scatterings are described by a single 
potential strength V. For a low impurity concentration rij, the quasiparticle self-energy 
in the SCTMA |2^] can be expressed in the Nambu spinor representation as 

T 1 ^ A \e)=n l T R{ - A \e) = (A^^o+TO, (2.1) 

for |e| <C 7- Here A is the mass renormalization factor, 7 is the impurity-induced relax- 
ation rate, r/ is a dimensionless parameter, 777 represents the decrease of the chemical 
potential induced by the impurity scatterings, and tq and Tj (i = 1, 2, 3) stand for the 



2x2 unity and Pauli matrices, respectively. A use of Dyson's equation immediately 
yields the impurity-averaged one-particle Green's functions as 

k () " [(l-A)6±i 7 ]2- £ 2 ' 

where the shift of chemical potential has been absorbed in jjl. The zero-energy density of 
states is calculated as po = — (l/7r)Im ^ fc T?G R = Alj/ir 2 VfV g , where G R ^ = G R ^ A \o) 
and I = ln(r/7) > 1 with F ~ ^jVfVg/a. 

The parameters 7, A, and 77 can be evaluated consistently via the T-matrix equation 

T Ri - A \e) = Vt 3 + Vr 3 g R{A) {e)T R ( A) (e), (2.3) 

with g R{A \e) = £ fc Gj (A) (e). Using Eq. (2.2), we can show that 

g R(A) {e) = m [(A - 1)(1 - r*)e T 17] + (U- 1 - CT 1 ) ts, (2.4) 

for |e| <C 7, where U is the effective impurity potential given by U~ l = V^ 1 + 
EfeW(4 + 7 2 )- A substitution of Eqs. (2.1) and (2.4) into Eq. (2.3) leads to 7 = 
2^/71700(1 + V 2 ), A = (1 - r, 2 ){l - 1)/( V 2 + 21- 1), and v = 2/irp U. 
The Born limit corresponds to rj 2 3> 21, yielding 

7 = ^n iPo U 2 , A = 1 - /; (2.5) 

and the unitary limit corresponds to rj — > 0, meaning that 

7 = — > A = - — -. (2.6) 
7rpo zt — 1 

Throughout this paper, we mainly consider the cases near either the Born or the unitary 
limit. It is worthy to point out that the values of the impurity potential V driving the 
system into these two limits are dependent on the band structure. 



III. The diffusive modes 

Since the QI effect is related to the diffusive modes, we first derive the expressions 
of the 0-mode and 7r-mode cooperons and diffusons for a disordered <i-wave supercon- 
ductor. 
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A. O-mode cooperon and diffuson 



The O-mode cooperon and diffuson exist both in RA and in RR channels due to 
the local particle-hole symmetry T2G R (e)T2 = —G^(—e). The ladder diagrams for the 
cooperon are given by Fig. 1(b) in Ref. ^3|. The equation for O-mode cooperon can be 
expressed as 

C(q; e, e') = W(e, e') + W(e, e')H(q; e, e')C(q; e, e'), (3.1) 

where the two-particle irreducible vertex W(e, e') and the integral kernel TC(q; e, e') are 
defined in the RR and RA channels as 

W(e, e ') RR ^ = ni T R (e) ® T R ^ (e ; ) (3.2) 

and 

H{q;e,e') RR{A) = ^G R _ k (e) <g> G R{A) (e). (3.3) 

k 

Equation (3.1) can be also expressed as 

A(q;e,e')C(q;e,e')=W(e,e'), (3.4) 

where 

A(q; e, e')=l- W(e, e')H(q; e, e'), (3.5) 

with 2 = To (g> To. From Eqs. (3.1)-(3.5), it follows that one can make a decomposition 
of X = Y.ij X ij T i ® T j for x = W,H,A, and C. Substituting Eq. (2.1) into Eq. (3.2), 
we can easily obtain all the nonvanishing components of W(e, e') as 



wp (l +rf) 



l + i~(e±e') 

7 



(3.6) 
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and 



l±i-e' 



(3.9) 



;3 ° ~ T \po(l + ^)\ 7 
It then follows that the dominant component of W(e, e') near the Born limit (r/ 2 » 21) 

is W(e,e')s3, while that near the unitary limit (rj 2 <C 1) is W(e,e')oo- 

Now let us calculate H(q;e,e') and *4(g;e, e'). For small values of q, e, and e', we 

have 



Gj_ fc (e) « G<£ + e^G H (e') | e , =0 -g • VG<£ + ^gg : VVG£ (3.10) 
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1 



and 



de 



G fc l e ) ~ G fc + e ^ G fc l e J U=o • 



(3.11) 



Substituting Eqs. (3.10) and (3.11) into Eq. (3.3), and using Eqs. (A1)-(A4) in the 
Appendix, we can show that all the nonvanishing components of H(q;e,e') are given 
by 



4/ 7 



,2 Uy 2 



and 



On -g^ 



H7„- * e '\ RR ( A ) - n P° 



(2Z-l)+i^-^(e±e / ) 1 



7 



12 7 2 



(2/ - 1) +i (e±e') - 



<7 



(3.12) 
(3.13) 

(3.14) 



7 ■ 12 7 2 

The diagonal components of A(q; e, e') can be easily calculated by a substitution of 
Eqs. (3.6), (3.7), and (3.12)-(3.14) into Eq. (3.5), yielding 

1 



A(q;e,e'C (A) = 



4Z(1 + ?7 2 ) 



2(2/ - 1) - (21 + l)r] 2 



+i 



,2A + (1-A)r/ 2 



(e ± e>) - 



2v 



2 + 2v 2 + (3vj + v 2 ) n 
12 7 2 



(3.15) 



A(q;e,e>)« R W=± 



4/(1 + V 2 ) 



.2A(Z-1) + 1^^ 

v2 * 



(2i-l)+i=^! ^_J_l( e ± e ') ■" ■ 
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12 7 2 



,(3.16) 



A(q;e,e>)« R ^ - > 



4/(1 +7/ 2 ) 



1 - A , ^ ? + 3f 2 
( 2/-l ) + i^ r ( e ±e0-^^g 2 



(3.17) 
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and 



A(q;e,e') RR ^=±-^[(2l-l) + 2r ] 



+i 2A "-" + 1 (e ± /) - t » +3 "' + 1 ff + 2 "» ) " 2 ^ . (3.18) 



We note that all the non-diagonal components of H(q; e, e') are vanishing, and those 
of W(e, e') are negligible near the Born or unitary limit. Then Eq. (3.1) or (3.4) indicates 
that only the diagonal components of C(q; e, e') may be singular for these two limits. As 
a result, the 0-mode cooperon can be expressed as C(q; e, e') = J^i C(q; e, e')ii r i ® T i- As 
will be shown below, all C(q; e, e')u are of diffusive poles. Obviously, only the diagonal 
components of A(q; e, e') and W(e, e') are needed for the calculation of 0-mode cooperon, 
and thus Eq. (3.4) becomes equivalent to the following group of equations: 

AxiCoo + AuCu + A 22 C 2 2 + A 33 C 33 = Wqo, (3.19) 

A 00 C U + ^nCoo - A 22 C 33 - A 33 C 22 = 0, (3.20) 
A)C)C22 — ^4llC33 + A22C00 — A 33 C\\ = 0, (3-21) 
A)()C33 — A11C22 — A22C11 + ^33600 = W 33 , (3.22) 

where the arguments (q, e, and e') of An, C„, and Wu have been omitted. 

Near the Born limit (rj 2 3> 21), the singular terms in Eqs. (3.19)-(3.22) satisfy the 
following relations: 

(21 + 1)C RR(A) + (21 - 1)C RR(A) ± 2C RR[A) = 0, (3.23) 

(21 + 1)C RR{A) - (21 - l)C RR{A) T 2C RR{A) = 0, (3.24) 
(21 + l)C RR(A) + (21 - 1)C RR(A) T 2C RR(A) = 0, (3.25) 
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(21 + 1)C RR{A) - (21 - 1)C™ (A) ± 2C fi fl(A) = 0, (3.26) 

the only solution of which is given by 

n RR{A) _ n RR{A) _ n RR{A) _ ^ n RR(A) ,„ 
°00 r°ll — _ °22 ^33 • 

Substituting Eqs. (3.7), (3.15)-(3.18), and (3.27) into Eq. (3.22), and using Eq. (2.5), 
we can show that the terms of order r]~ 2 cancel exactly out, leading to 

C(q;e,e')« RiA) = T^ n2 \ , , y (3.28) 
7rpo Dq z — i(e ± e) 

where D = (v 2 + v 2 ) / Alj is the quasiparticle diffusion coefficient. Combining Eq. (3.27) 
with Eq. (3.28), and noting that the 0-mode diffuson has the same expression as that 
of the 0-mode cooperon due to the time-reversal symmetry, we obtain 

C(q-e,d) RR W =V(q;e,e') RR W = ^- 1 , J tto®t +ti®ti±7s®75+T3®t 3 Y 

7rpo Dq — t\£ ± e'J \ / 

(3.29) 



which is in agreement with that appearing in Ref. |2^] (aside from a disputed pre- 
factor). Similarly, we can show that Eq. (3.29) is also valid near the unitary limit 
(i] 2 <C 1). The above evaluations indicate that any small deviations from either limit 
do not make the Goldstone 0-modes gapped. However, in the intermediate region far 
from these two limits, the non-diagonal components of W(e, e') cannot be neglected, 
and thus C(q;e,e') may contain some singular non-diagonal components. 

B. 7r-mode cooperon and diffuson 

In the UN limit, there exist the additional 7r-mode cooperon and diffuson due 
to the global particle-hole symmetry |2^, 26] T2G R ^ A \e)T2 = G^^^e), with Q = 



±(7r/a, ±7r/a) the nesting vector. Any small deviations from this combined limit can 
be shown to make the Goldstone w- modes gapped. The ladder diagrams for the 7r-mode 
cooperon can be obtained from those of the 0-mode cooperon by replacing q by Q + q. 
The equation for the 7r-mode cooperon is given by 

A w (q; e, e')C n (q; e, e') = W(e, e'), (3.30) 



where 

Arte; e, e') = X - W(e, e')^(q; e, e'), (3.31) 

with 

H„(q;e,e') RR{A) = ^2 GQ +q _ k {e) ® G^ {A \e). (3.32) 

k 

In order to calculate 7i n (q; e, e') and A n (q; e, e'), one needs to exploit the relations 
^Q+fc = and £,Q + k = —£,k — 2/i. For a nearly-nested Fermi surface (|/i| <C 7), we 
have 

C R (A ^ tC r (At , o f 2C fc (»7ro - A fc n) + ( 7 2 + A| - ^) T3 



+ (»)' h ' + A * - ^ )( '^ T A "> - ^ + 3A ^ - (i)T3 . (3.33) 

17 + e fc) 

Substituting Eqs. (3.11) and (3.33) into Eq. (3.32), and using Eqs. (A1)-(A4), we can 
show that all the nonvanishing components of H w (q; e, e') are given by 

TT In- , J\ RRi - A ) - U(„. , J\ RRi - A ) + 71 ' P °^ (<\ <\A\ 

^7T{q,e,e) 00 -H(q,e,e) 00 ± 3 , (6.64) 

tt i i\RR(A) tji i\RR(A) . npo/J 2 , s 

H w (q; e, e ) n = -H(q; e, e ) n + 7^3-, (3.35) 

„ / r^RR(A) tji i\RR{A) . npo/j 2 

Hn(q; e, e ) 3 3 = ~H(q; e, e ) 3 3 + ' ( 3 - 36 ) 



and 



„ i f\RR(A) .tt, ,^RR(A) .7T/?oM /9 q m 

^(91 e, e ) 03 = ±H n {q; e, e ) 30 = (3-37) 



A substitution of Eqs. (3.6)-(3.9) and (3.34)-(3.37) into Eq. (3.31) yields the diagonal 
components of A w (q; e, e') as 

- /\jr 2 

) 

,2 , o„,2 fq„,2 , „,2\„2 



AM: e, e'C (A) = "iTTT^y [ - 2(2/ - 1) - (6/ - lW + i 2A ~ ~ A ^ ( e ± e ') 



2^ + 2^ - (3^ + ijV 2 (4 - 6 V 2 )fi 2 4 W 
12 7 2 9 3 7 2 7 



}3.38) 



2 ! . S „2 „2. 



A u (g; e, e ) n = T ^ + ^ [(21 - 1) + 1 (e ± e ) - ^ q - ^ {3.39) 



7 v ' 12 7 2 * 37 2 
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A n (q;e,e) 22 ' = - — - [(2/ - 1) + ,— (e ± e ) - -J^Q - , (3-40) 
and 



41(1 + rj 2 ) 

vj + 3^ - (2vj + 2t; 2 )r/ 2 2 4 W 



1 

\ / / 2 Ann n~\ 

(3.41) 



12 7 2 " 37 2 7 

Like the 0-mode cooperon, the 7r-mode cooperon near the Born or unitary limit can 
be also expressed as C n (q; e, e') = Yli Ctt(q; Oi* r » ® r i> with all C n (q; e, e')^ assumed 
to be of diffusive poles. Therefore, Eqs. (3.19)-(3.22) are also suitable for the 7r-mode 
cooperon. Near the unitary limit (rf <C 1), the singular terms in Eqs. (3.19)-(3.22) 
(for the 7r-mode cooperon) satisfy the following relations: 

ori RR{A) s~,RR(A) ri RR{A) _ „ /„ A0 \ 

ZU tt00 t Wll T W33 — u > l"3- 4Z J 



2C RR(A) T C ^(A) ± ^(A) = Qj (3>43) 



2Cf 2 ? A) ± C*S A) ± eg? A) = 0, (3.44) 



on RR(A) , r ,RR(A) n RR.(A) _ „ ,„ . _s 

ZU tt33 ± °7r22 T U tt00 — U ' ^.4DJ 

the solution of which is shown to be 

n RR{A) _ , n RR(A) _ n RR{A) _ , n RR(A) ,„ 

°7r00 - ±0 7rll - _C/ 7r22 ~ ±L/ 7r33 ■ (3.4b j 

Substituting Eqs. (3.6), (3.38)-(3.41) and (3.46) into Eq. (3.19), and using Eq. (2.6), 
we can show that 



j^RR(A) _ 1_ 

ixpoDq 2 -i(e±e') + 25' 



C^q; e, e')oo 1 j = T-r n , ~ , ^ , OJC , (3.47) 



with 5 = 2r/ 2 7 + 277^/7 + /j, 2 /Ij <C 7. The present expression of (5 supports the result 
estimated by Yashenkin et al. [^3|. Combining Eq. (3.46) with Eq. (3.47), and noting 
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that the 7r-mode diffuson has also the same expression as that of the 7r-mode cooperon 
due to the time-reversal symmetry, we obtain 

C^q;e,e') RR{A) =^e,e') RR{A) 
47 2 1 / \ 

= 771 ., T A , n A ± r ° ® r » + r l ® r i =F T 2 ® T 2 + T 3 ® ^"3 • (3.48) 

Near the Born limit (r/ 2 3> 21), one can easily show that Cf^" = 0, (i = 0,1,2,3), 
indicating that the diffusive 7r-modes exist only near the UN limit. Contrary to the 
diffusive 0-modes, the Goldstone 7r-modes are gapped by any small deviations from the 
UN limit measured by 5. For the situations far from the UN limit, the contributions of 
diffusive 7r-modes to the QI effect are completely suppressed due to the large gap. 

IV. QI correction to the quasiparticle DOS 

The QI correction to the quasiparticle DOS can be calculated via 

A/0 ( e ) = -Ilm^TrAG^(6), (4.1) 

k 

where AG R (e) represents the lowest-order correction to the one-particle Green's func- 
tion due to the diffusive modes. 

A. Contribution of the 0-mode cooperon 

The lowest-order self-energy diagrams with 0-mode cooperon are depicted in Fig. 1. 
The contribution of Fig. 1(a) has been calculated in Ref. [p3; ]. Figure 1(b) contains the 
non-singular ladders, its contribution to the DOS can be shown to be the same order 
of magnitude as that of Fig. 1(a). According to Eq. (4.1), the contribution of Fig. 1(a) 
to the DOS reads 

p(e) a = -ilm^^C(q;e,e)^Tr (nGq_ k nG^G R ) . 

kq i 
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(a) (b) 

Figure 1: The lowest-order self-energy diagrams with 0-mode cooperon (shaded blocks). 
The grey block in Fig. 1(b) represents the non-singular ladders. The corresponding 
diagrams with 7r-mode cooperon can be obtained by the replacement of q — + Q + q in 
Figs. 1(a) and 1(b). 

By means of Eq. (A5), it can be rewritten as 

p(e) a = -^Im^^[C( g ;e, e )^ a ] i ., 

q i 

with M a = G R ®(G R G R ) . The contribution of Fig. 1(b) can be similarly evaluated, 
and the total contribution of Figs. 1(a) and 1(b) can be expressed by 

p(6)coo P = ~Im EE [C(Q;e,e) RR M] u , (4.2) 
q i 

with 

M = Y,G R ®{G R KG R ). (4.3) 
k 

Here the matrix K satisfies 

K = t + n % E T R G R KG R T R , (4.4) 
k 

with T R ^ = T R( - A \0). Substituting Eqs. (2.1) and (2.2) into Eq. (4.4), and using 
Eqs. (A1)-(A4), we can show that 

K = (VC)ro, (4.5) 
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where 

1, for the Born limit, 



c 



(4.6) 

21 — 1, for the unitary limit. 



Equation (4.5) leads to M. = (l/QM. a , and hence p(e) coop = (l/()p(e) a , indicating that 
the contribution of Fig. 1(b) just renormalizes the pre-factor of p(e) a . A substitution 
of Eqs. (2.2) and (4.5) into Eq. (4.3) yields 

M = -i^—r 1 (n <g> n + r 3 <g> t 3 ) . (4.7) 



27rCu/v 5 7 

Substituting Eqs. (3.29) and (4.7) into Eq. (4.2), and noting that the upper cutoff of q 



can be set to be l/l e with l e = D /2~f the mean- free path, we obtain 

p(e)coo P _ 27r T3„Y^ D 



po a( ^ Dq 2 - i2e 



1 Inrr ( 4 - 8 ) 



2a( \e 
where a = (v 2 + v 2 )/2vfV g . 

B. Vanishing contribution of the O-mode diffuson 



The lowest-order self-energy diagrams with the O-mode diffuson are given by Figs. 2(a) 
and 2(b), and the vertex correction to the impurity-scattering T-matrix is shown by 
Fig. 2(c). Unlike the single- vertex-dressed diagram used in Ref. J23|, the present 
Figs. 2(a) and 2(b) include the additional diagrams with both vertices dressed by the 
7r-mode diffuson. Similar vertex correction has been considered in the theory of dis- 



ordered interacting-electron systems 2£]. The vertex-corrected retarded T-matrix can 
be expressed by 



f R (q, e)^ = ^ J(q, e)™^ u ,T R (e) uu ,, (U 



VI' 



where fi,fx',u, and v' are the indices of the Nambu space, and the vertex function 
J{q, e) RR is given by 

J(q, e) RR =1 + V(q; e, e) RR H(q; e, e) RR . (4.10) 
14 





k' k q+k k W 






(c) 



Figure 2: The lowest-order self-energy diagrams with 0-mode diffuson (a and b), and 
the vertex diagram (c). The shaded and grey blocks represent, respectively, the 0-mode 
diffuson and non-singular ladders. The dashed lines denote the impurity-scattering T- 
matrix. The corresponding diagrams with 7r-mode diffuson can be generated by the 
replacement of q — + Q + q in these 0-mode diagrams. 



In order to calculate the vertex function, we exploit the equation for 0-mode diffuson 
in the RR-channel 



V(q; e, e) RR = W(e, e) RH + W(e, e) HH H{q; e, e) KH V(q; e, e) HH , (4.11 



\RR 



\RR 



RR< 



\RR 



yielding 



H{q; e, e) RR = W~\e, e) RR - V~ l (q- e, e) RR . 
A substitution of Eq. (4.12) into Eq. (4.10) leads to 

J(q, e) RR = V{q; e, e) RR W-\e, e) RR . 

Making use of Eqs. (3.6) and (3.7), we get 



2 7 



r 3 ® T3, for the Born limit, 
-tq tq, for the unitary limit. 



(4.12) 



(4.13) 



(4.14) 
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Substituting Eqs. (3.29) and (4.14) into Eq. (4.13), we obtain the expression of the 
vertex function as (for Dq 2 <C 7 and |e| <C 7) 

J(q, e) RR = 2 ^ — (r ® t - n <g> ti - t 2 ® r 2 - r 3 <g> r 3 ) , (4.15) 

which is suitable both near the Born and near the unitary limits. Substituting Eqs. (2.1) 
and (4.15) into Eq. (4.9), we can easily show that 

f R (q, e) = £ J(q, e^T^r* = 0, (4.16) 

i 

indicating that both Figs. 2(a) and 2(b) have vanishing contributions to the DOS near 
the Born or unitary limit. 

C. Contribution of the 7r-mode cooperon near the UN limit 

Near the UN limit, besides the diffusive 0- modes, the Goldstone it- modes also con- 
tribute to the QI effect. The leading self-energy diagrams with the 7r-mode cooperon 
can be obtained from those in Fig. 1 by a replacement of q — ► q + Q. Similarly with 
the case of 0-mode cooperon, the contribution of 7r-mode cooperon to the DOS can be 
evaluated via 

q i 

where 

M n = Y,G^. k ^(G R KG R ), (4.18) 
k 

with K = Itq/(21 — 1) for the unitary limit. Expression (4.18) can be readily evaluated 
by means of the global particle-hole symmetry, yielding 

M « = 1 1 2,(21 -1) W ^ Tl + T3 ^ T3 ) • (419) 
Substituting Eqs. (3.48) and (4.19) into Eq. (4.17), we obtain 
/o(e) 7r _ COO p _ 2i ^ D 



M2T^T) Re ^ 



po a(2l - 1) ^ Dq 2 - i2e + 25 



In JL- . (4.20) 



2a(2£ - 1) ^TS 2 " 
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Comparing Eq. (4.8) with Eq. (4.20), one finds that the contribution to DOS of the 
7r-mode cooperon has an equal magnitude but an opposite sign to that of the 0-mode 
cooperon in the UN limit (<5 —►()). 

D. Contribution of the 7r-mode diffuson near the UN limt 

The leading self-energy diagrams containing the 7r-mode diffuson are generated by 
replacing q by Q + q in Figs. 2(a) and 2(b). Obviously, Eqs. (4.9) and (4.13) are also 
suitable for the 7r-mode diffuson. By exploiting Eq. (3.48), we can readily show that 

,RR 2 7 



and 



Jw(q, e) = .' , ns T <g> T + n ® ti - t 2 <8> r 2 + t 3 <g> r 3 (4.21) 

— z2e + 2d 



% 

I67 1 , , 

= -^ ZV-*2 £ + 2^ (422) 

which are valid near the UN limit. 

The contribution of the 7r-mode diffuson to the DOS is given by 



p(e).-difr = Tr z) G Q+kT?(q, e)G*KG%] . (4.23) 

kq 



Substituting Eqs. (2.2) and (4.22) into Eq. (4.23), and using Eqs. (2.6) and (Al)-(A4), 
one can easily show that 

l(e)7r-diff 32tt7 \ - D 

p a(2Z-l) e ^ (Dq 2 - i2e + 25) 2 

4 7(5 (4.24) 



a(2Z-l)e 2 + J 2 ' 



Expression (4.24) is quite different from the logarithmic behavior obtained in Ref. 23] 
due to the additional contributions of the self-energy diagrams with both the impurity 
scattering vertices dressed by the 7r-mode cooperon. 
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E. Results of the QI correction to DOS 



For generic Fermi surfaces, the QI correction to the DOS results only from the 
contribution of 0-mode cooperon, i.e., 

Mf) = ^)coo P= _J_ ln 7 (425) 

Po Po 2aC |e| 

where £ is given by Eq. (4.6). Therefore, the QI effect gives rise to a logarithmic 
suppression to the quasiparticle DOS, as predicted by Yashenkin et al. Equa- 
tion (4.25) is suitable near the Born or unitary limit. In the intermediate region far 
from these two limits, however, the 0-mode cooperon may contain some singular non- 
diagonal components, and thus a refined theory is necessary for the generic situations. 
In addition, the additional contributions of the self-energy diagrams with non-singular 
ladders yield various renormalization factors in these two limits. 

Near the UN limit, the QI correction to the DOS is the sum of p(e) COO p, /o(e) 7r _ COO p, 
and /9(e) 7r _difr, so that 



Ap(e)uN 1 



po 2a(2l - 1) 



lnP; + ln^= + 875 



(4.26) 



V^TP e 2 + 5\ 

The above equation indicates that the quasiparticle DOS near the UN limit can be 
subject to a positive correction due to the contributions of the Goldstone 7r-modes. 
Upon approaching the UN limit (6 is small enough), the contributions from the 0-mode 
and 7r-mode cooperons cancel out each other, so that the DOS correction is given only 
by the contribution of the 7r-mode diffuson, i.e., 

lim Ap( C ) UN = 47T7 () (427) 
This result can be used to account for the sharp peak found in the numerical studies 



of the DOS [19,2C|. 



V. Summary 
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We have extensively investigated the QI effect on the quasiparticle DOS in a d x i_ y 2- 
wave superconductor with dilute nonmagnetic impurities. As in the study of disordered 
normal metals fl24j| , the understanding of the diffusive modes is essential for the inves- 
tigation of the QI effect in a superconductor. Through detailed derivations, we have 
obtained the expressions for the diffusive modes both near the Born and near the uni- 
tary limits. It is demonstrated that these Goldstone modes may contain non-diagonal 
components in the intermediate region far from these two limits. Therefore, a refined 
weak- localization theory is necessary for the generic situations. 

For generic Fermi surfaces, the QI effect results only from the 0-mode cooperon, 
yielding a logarithmic suppression to the quasiparticle DOS near the Born or unitary 
limit. We show that those non-trivial self-energy diagrams containing non-singular 
ladders give rise to various renormalization factors of the DOS correction for the Born 
and unitary limits. 

Near the UN (unitary and nesting) limit, the QI effect comes not only from the 
contribution of the 0-mode cooperon, but also from those of the 7r-mode cooperon 
and diffuson. It is found that the self-energy diagrams with both impurity-scattering 
vertices corrected by the 7r-mode diffuson have an important contribution to the DOS, 
which is proportional to 5/(e 2 + 5 2 ). As a result, the quasiparticle DOS is subject to a 
positive correction induced by the diffusive 7r-modes. Upon approaching the UN limit 
(<5 — ► 0), the contributions of the 0-mode and 7r-mode cooperons cancel out each other 
(even including the renormalization contributions induced by the non-singular ladders) , 
so that the DOS correction becomes a 5-function of the energy. This result can be used 
to account for the resonant peak found in the previous numerical studies |l^, |20j. A 
similar sharp peak of electronic DOS has been found numerically in a disordered 2D 
tight-binding model for the normal state |3(|. The appearance of this sharp peak can 
be also explained by the contribution of Fig. 2(a) with the 7r-mode diffuson (Fig. 2(b) 
has a vanishing contribution for the normal state) . 

Like in the case of disordered interacting-electron systems [^| , all the leading polar- 
ization diagrams responsible for the QI effect on the quasiparticle conductitity can be 
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generated in the conserving approximation from the lowest-order self-energy diagrams 
shown in Figs. 1 and 2. For the normal state, it has been shown that the contributions 
of those polarization diagrams with 7r-mode diffuson lead to an antilocalization correc- 
tion to the conductivity in the UN limit [32]. How the QI processes related with the 



diffusive 7r-modes affect the transport properties, such as electrical and spin conductiv- 
ities, in a d-wave superconductor is another interesting problem, and will be studied in 
our future work. 
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Appendix: Some useful mathematical formulas 

In this appendix, we present some mathematical formulas, which are useful for the 
evaluations in the previous sections. By the approach used in Ref. [33|, we can show 
that 

F 2 /\ 2 1 
^( 7 2 + e fc) 4 2 4nv f v g ^ K ) 

V & _V" A fc 1 ( A2) 

V^ 2 + 4) 4 V^ 2 + 4) 4 8tt W 2 ' 1 ' 

?(7 2 + 4) n = (n-l)^ W 2(«-D' forn ^ 2 ' ( A3 ) 



and 



V ^fc - V ^fc t for n > 3 (Ad) 



20 



As an example, we shall prove Eq. (Al). Using the Dirac-type quasiparticle spectrum, 
and noting that there exist four gap nodes, we have 

£ A S _ A f f d kdk vjvfkfkl 



y UK A [ fdkfc 

V (r 2 + 4) 4 J J 



By means of the transformations of pj = \/vj /v g kf and p g = y/vg/v/kg, the above 
equation can be changed as 

V + i) 4 ■/ ■/ ( 27r ) 2 fr 2 + »/ w «(pJ +^)] 4 

= — 7. n I dx- -7 / d6 cos 2 #sin 2 9, 

2^v f v gl i J (1 + xY Jo 

where xq = VfVgp^j^ 1 with po ~ I /a. For the weak-disorder case (7 is small enough) 
under considered, we can set xq = 00. Thus, a completion of the integrals over x and 
9 in the above equation immediately yields Eq. (Al). 
Another useful formula is given by 

X - CuTrinAnB) = Y^(CM)a, (A5) 

i i 

where C = J2i Cun^Ti, and M = A®B with A and B being any linear superimpositions 
of Ti (i=0, 1,2,3). In order to prove Eq. (A5), we assum that T{A = J2j x ij T j an d 
nB = VikTk- Then we have 

- ^2 CuTiinATiB) = 2^2 CiiXijyikTviTjTk) = ^ (',,■>■, jH.i.^ji, = ^ <"„•<•,., .</,., (A6) 

i ijk ijk ij 

and 

CM = Y^ CuinA) ® ( n B) = ('„■>■; , mi,- ® r fc . (A7) 

i ijk 

From Eq. (A7), it follows that 

^2{CM)u = Y CnXijVij. (A8) 

i ij 

A combination of Eq. (A6) with Eq. (A8) immediately leads to Eq. (A5). 
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